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Abstract 

In this report, we present a systematic account of mathematical structures of certain special 
polynomials arisen from the energy study of the superintegrable TV-state chiral Potts model with 
a finite number of sizes. The polynomials of low-lying sectors are represented in two different 
forms, one of which is directly related to the energy description of superintegrable chiral Potts 
ZAT-spin chain via the representation theory of Onsager's algebra. Both two types of polynomials 
satisfy some (A'^-l-l)-term recurrence relations, and iVth order differential equations; polynomials 
of one kind reveal certain Chebyshev-like properties. Here we provide a rigorous mathematical 
argument for cases A'^ = 2, 3, and further raise some mathematical conjectures on those special 
polynomials for a general N . 
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1 Introduction 



This article contains a sound mathematical treatment of certain special polynomials appeared in the 
eigenvalue problem of superintegrable chiral Potts Z^v-spin chain of a finite size through Onsager's 

27ri 

algebra representation. Throughout this paper, N will always denote an integer > 2, a; = e jv , and 

is the vector space consisting of A^-cyclic vectors with the basis {|'^)}„gZjv' Denote by X,Z 

the operators of satisfying the Weyl relation and A^th power identity property: XZ = u}~^ZX, 

X^ = Z^ = 1. One can represent X,Z in matrix forms by X\n) = |n + 1), Z\n) = oj'^\n) for 

n G Zjv- The Hamiltonian of the superintegrable chiral Potts Zjv-spin chain is the operator acting 

L ^ 
on (g) C defined by 

Hik') = - E E + k'zrz[l-n (1) 

1=1 n=l 

where k' is a real parameter, and Xi, Zi are operators X, Z acting on at the site I with periodic 
boundary condition: Zl+i = Z\. For N = 2, this is the Ising quantum chain. The Hamiltonian 
H{k') commutes with the spin-shift operator ^i, hence carries the ZAr-charge, denoted hy Q = 
0, . . . , A?" — 1, in this paper. One can split H{k') into a sum of Dolan-Grady pair of operators with 
the coupling parameter fc', by which the representation theory of Onsager's algebra can be employed 
for the study of -eigenvalue problem [13]. However for its spectrum calculation, the available 
method relies on the iV-state chiral Potts model, a two-dimensional solvable statistical lattice model, 
where H{k') is derived as the derivative of transfer matrix at the superintegrable rapidity point 
(see, e.g. [1, 3] and references therein). For the superintegrable chiral Potts model with a finite size, 
the eigenvalues of transfer matrices are given by certain polynomial solutions of a spin-^^^^^ XXZ- 
model-like Bethe equation [1, 3, 4, 5, 7]. The low-lying (no Bethe-excitation) sector polynomials 
of chiral Potts transfer matrices for the ground-state eigenvalue were explicitly constructed in [1, 
3, 5], and will be called the Baxter- Albertini-McCoy- Perk (BAMP) polynomials in this article. By 
converting the chiral Potts eigen-polynomials to another ones through a simple change of variables, 
one can determine the eigenvalues of the Hamiltonian (1). In particular, the BAMP polynomials 
are transfered to another type of polynomials, which are more convenient to express the spectra 
of (1). The resulting polynomials in the Ising case turn out to be Chebyshev polynomials. For 
> 3, these polynomials have shown many properties reminiscent of, but not fully in agreement 
with orthogonal polynomials, hence will be called the (generalized) Chebyshev-type polynomials 
in this paper. In [12, 14], we have made empirical calculations on these polynomials for small N 
by investigating the recurrence relations, differential equations etc., in hope to explore the zero 
distributions of those Chebyshev-type polynomials. The computations have shown a remarkable 
feature of "partial orthogonality" in terms of Jacobi weights, leading to a certain conjecture about 
those Chebyshev-type polynomials (for the details, see [12]). Through explicit calculations made in 
the study of these polynomials in the context of superintegrable chiral Potts models, we discovered 
certain simple mathematical structures in the process. It is the aim of this report to present a 
mathematical account on the structures of the BAMP and Chebyshev-type polynomials for all 
N regarding to the recurrence relations and differential equations, based upon the computational 
results found in [12, 14]. The clarified mathematical presentation of the subject in this article, we 
hope, could stimulate some interest in these special polynomials with a sound physical interest, 
however unexplored in classical literatures ( e.g. [6, 17]) to the best of our knowledge. 

The paper is organized as follows. In Sec. 2, we review some basic facts on Onsager's algebra 
and its connection with the Hamiltonian (1). In Sec. 3, we define the BAMP polynomials appeared 
in the spectra of A^^-state chiral Potts models, and discuss their common mathematical structures 
for all A''. In the Ising case, the relation between ground state energy and all other eigenvalues are 
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given. For a general N, the recurrence relation of BAMP polynomials by varying the size L modular 
A'^ is systematically constructed, and certain qualitative properties related to the recurrence relation 
are discussed. The problem of hypergeometry-like higher order differential equations governing the 
BAMP polynomials are raised in Sec. 4 with illustrative discussions on the case N = 2,3. In 
Sect. 5, we define the Chebyshev-type polynomials, and discuss their recurrence relations. In the 
Ising case, we provide a detailed account on its relation with Chebyshev polynomials. In Sect. 
6, we discuss the problem of differential equations for Chebyshev-type polynomials, and present a 
rigorous mathematical proof of solutions for those "Z3- Chebyshev- type" polynomials. 

2 Onsager's Algebra 

In the seminal paper of Onsager on the free energy solution of the two-dimensional Ising model 
[15], there appeared an infinite-dimensional Lie algebra, called Onsager's algebra, generated by a 
basis {Am,Gi}jn^z leN ^^^^ the following commutation relations: 

[■^m, ^n] = ^Gm-n , [-^m, Gi] = 2Am_i — 2Am^i , [Gm, Gi] = , 

where G-i := —Gi and Go := 0. The elements Aq,Ai satisfy the Dolan-Grady relation [11]: 

[Ai,[Ai,[Ai,Ao]]] = 16[Ai,Ao] , [Ao,[Ao,[Ao, A,]]] = 16[^o,^i] , 

and Onsager's algebra is the Lie algebra generated by this Dolan-Grady pair {^0)^i} [10, 16]. 
Inspired by results in [9] on the (finite-dimensional) unitary representations of Onsager's algebra, 
a realization of Onsager's algebra as the Lie-subalgebra of the loop algebra sl2[z,z~^] fixed by a 
standard involution was found in [16] through the identification: 

Am = 2z"'e+ + 2z-"'e- , Gm = {z"" - z-'^)h , m G Z, 

where e^,h are generators of sl2 with [e+,e~] = h, [h,e^] = ib2e^. By which, a thorough math- 
ematical study of Onsager's algebra was carried out in [8] by using techniques in ideal theory 
of polynomial algebras, and the discussions further enrich our understanding the various aspects 
of the mathematical structures of Onsager's algebra. It is known that all finite-dimensional ir- 
reducible representations of the loop algebra b12[z,z~^] are given by tensoring a finite number 
of irreducible si2-representations through the evaluation of z at distinct non-zero complex values 
ajS. The Hermitian irreducible representations of Onsager's algebra are obtained by factoring On- 
sager's algebra through the s/2[-2, ^~^]-representations, with further constraints on the evaluated 
values ttj's: \a.j\ = l,aj ^ ±1 and aj / a^^ for j / k. By which. Am = 2X]j=i(a™e^ + '^J"^^j)i 
Gm = Z]j=i(flj" ~ aj"'')hj. For the Hamiltonian (1), we write H{k') as a sum of Dolan-Grady pair 
{Ao;^i} with the parameter k': H{k') = — ^(^o + k'Ai). By [3], all the sZ2-irreducible represen- 
tations involved in the spectra of H{k') are of spin-^, hence their eigenvalues have the Ising-like 
expression: 

n 

a + pk' + N ^±Jl-2k' cos{9j) + k'"^ , cos(^j) = ^{aj + aj^) , a, /3 G R . 

The number and location of Cj's, depending on the chain-size L and Z^v-charge Q, are the ingredients 
for the eigenvalues of H{k'). For the rest of this paper, we shall mainly discuss issues related to 
the ground-state eigenvalue. 
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3 Recurrence Relation of Baxter- Albertini-McCoy-Perk (B AMP) 
Polynomials 

We start with some simple elementary algebra for later use. For a given A^, the symbol m for m G Z 
will always denote the integer, < m < N — 1, with m = m (mod N). Let C[t], C[s] be the 
polynomial algebras in variables t, s respectively with the relation s := . Then any clement f(t) 
in C[t] can be uniquely expressed by f{t) = ^2^=0^ i'^ Pji^) with Fj{s) G C[s]. The multiplication 
by J^fSg^t^ on C[t], f{t) i-^ g{t) = {J2f=o^ i'')fi't)j gives rise to an endomorphism of C[t], which 
can be represented by the following matrix-form (by considering C[t] as a C[s]-module ) with 
9{t) = j:fJot^Gj{s), 



1 Go{s) \ 




f Fo(s) ] 




Gi{s) 








= R 


; R:= 


\ Gn-i{s) j 




[ Fn-i{s) ) 





1 1 



1 1 



V 1 



s 



(2) 



(3) 



Denote by A{x) the characteristic polynomial of the matrix : 

A{x) = det{x -R^)=x^ + ... + Ak{s)x^-'' + ... + An{s) . 

By the Hamilton-Cayley theorem, we have 

AiR^") = {R^'f + ... + Akis)iR^f-'' + ...+ Ajv(s) = . 

For an arbitrarily given initial polynomial b{t), we consider a sequence of s-polynomials Fij(s) 
for I G Z>o and < j < N — 1, defined by the relation: 



N-l N-1 
j=0 



(4) 



i=0 



Then A{R^) ^F; o(s), i(s), . . . , F; jv-i(s)^ = for Z G Z>o, equivalently, the following recurrence 
relation holds for s-polynomials Fij{s) with a fixed j, 



Fi+N,j{s) + ... + Ak{s)Fi+N-k,j{^) + ...+ An{s)Fij{s) = , l>0 . 



(5) 



Note that the same (A'^ + l)-recurrence relation for F^,j{s)^s in above holds for any given initial 
polynomial b{t). For N = 2,3, 4, the explicit forms of (5) are given by 



N 



F 



2{s + l)Fi+ij + {s- l)^Fij = 0, 



AT = 3, Fi+3,j - 3(s2 + 7s + l)Fi+2j + 3(s - l^Fi+ij - (s - ifFi,, = 0, 
AT = 4, F,+4j - 4(s + l)(s2 + 30s + 1)^^+3,,- + 2(35^ - 62s + 3)(s - lfFi+2,j 
-4(s + l)(s - + {s- 1)12^;,,- = 



(6) 



For a general AT, one can determine the values of A}^(1) for all k. Indeed, by setting s = 1 for the 

matrix R in (2), one can extend the vector (1, . . . , 1)* to a basis of so that Rs=i is expressed 
by dia[A'', 0, . . . , 0]. Hereafter, we use dia[aO) ■ ■ ■ > «Ar-i] to denote the diagonal matrix with entries 
aj. Hence the characteristic polynomial (3) of R^=i is equal to {x — N^)x^~^, which implies 



^i(l) 



-A^ 



N 



AkiX) = for k>2 . 



(7) 
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For the eigenvalue problem of the transfer matrix of size L in the superintegrable A^-state chiral 
Potts model, the eigenvalue of the ground state corresponds to the initial polynomial being 
the constant function 1. In general, the eigenvalues of the superintegrable chiral Potts model are 
solved by the method of Bethe equations in [5] as follows. Let /(i; v) be the i-polynomial of degree 
rup with the parameter v = {v\, . . . ^ v^p) defined by f{t; v) = Ytj^=i{l + ujvjt) . Associated to f{t; v), 
we consider the following function p{t; v): 

where Pa, Pb are integers satisfying < Pa + Pb < N — 1, Pi, — Pa = Q + L (mod N), with Q 
the Zjv-chargc as before. For an arbitrary v, p{t;v) is a rational function of t. The necessary and 
sufficient condition for p{t; v) to be a t-polynomial is the following Bethe constraints of -y, 

I -1 ™p -1 
\ _2 " "IT- j = l,...,mp. 9 

Vj +00 ^ Vj - UJVl 

Note that up to a phase factor the above equation is the Bethe equations for the spin- -^^^ , 
anisotropy 'J = jj XXZ chain with L-size and the periodic boundary condition [7]. For a Bethe- 
solution V of (9), p{t; v) is a i-polynomial, invariant under the transformation: 1 1— ut, hence it can 
be written as a polynomial of s (:= t^), 

p{t;v) = P{s;v) . 

Furthermore, the integers Pa,Pb are chosen so that ^(O;^) 7^ 0. The s-degree of the polynomial 
P{s;v) is given by deg P{s;v) = [ ^^ ^ ^"^'' ] with all its roots real, which determine the 
spectrum of (1). Here [r] denotes the integral part of a real number r. The eigenvalues of the 
ground state and one-excited state are the polynomials for rUp = and 1 respectively. 
For the Ising case, i.e. AT = 2, we have Pa + Pb = 0, 1. With f{t, v) = 1 in (8), one has 

= t') = ((1 + 0^ + - 0^) = ± n(- - s,) , 

whose zeros are known in [2]: sj = — tan^{J(j — ^— ^)}- By u) = —1, the right hand side 
of Bethe equations (9) is always equal to — Hence the solution of (9) is given by a 

collection of non-zero Vj's such that vj^ are roots of the above t-polynomial PjF!^p^. Hence the 
polynomial p{t,v) corresponding to a Bethe-solution of (9) has the form (up to a scale multiple): 

Pit, v) = Pis = t^v) ~ ±(n,,<...<i. s,) n;|,,,=i(^ - Si). 

For N > 3, the eigen-polynomial for the ground state of L-sizc transfer matrix in superintegrable 
chiral Potts model is determined by (8) by setting Pj, = 0, and fit;v) = 1 [1, 3], 

CW = Vg (l-t^t)^ ■ «,«:=zi^. (10) 

Note that aL,Q = 0'L',q ii L = L' (mod N). The Pq'\s) is a s-polynomial of degree 

.(iV-l)L-Q, 



bL,Q 



N 
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We shall call Pq'\s) the Baxter- Albertini-McCoy- Perk (BAMP) polynomial. For a given L, the 
polynomials Pq'\s), < Q < N — 1, are characterized by following relation: 



N-1 



(- 



= E t"''''PQ\^) ' ^L,Q = -L-Q (= (AT - 1)L - Q - NbL,Q) . (11) 
Q=0 



t - 1 



By which Pq'\s) are polynomials with positive integer coefficients. Furthermore, one has the 
following reciprocal property and recurrence relations among BAMP polynomials. 

Theorem 1 (i) Pq'\s) and Pq^\s) are reciprocal to each other if and only if L + Q + Q' = 
(mod N). 

(ii) Let Ak{s) he the polynomials defined by (3). "daa/tu d(^) 

recurrence relation by varying L for the same Q, 



The BAMP polynomials Pq satisfy the 



P^+^^)(.) + Ai(.)p(^+(^-^)^)(.) + . . . + ^,(.)pf +(^-'=)^)(.) + . . . + AMP^'\s) = . (12) 

(Hi) The polynomials Ak{s) are characterized by the property (12) for all L and Q. For 1 < k < 
N, Ak{s) is a reciprocal polynomial of degree k{N — 1) with the leading coefficient and Ak{0) equal 
to (-1)'=(^). Furthermore, Ai{s) = -NP^^\s), An{s) = (-l)^(s- 1)^(^-^, and for k> 2, one 
has Ay;{s) = {s — l)^''=^fe(s) with A^^l) ^ for some positive integer 1^. 

Proof By replacing t by | in (11), one obtains (^Jfr)^ = T,QZot^s''^''^P^^\l). Hence the 
relation P^\s) = s^i^'QP^\\) holds if and only if L + Q + Q' = (mod AT), then follows {i). 
For 0<m<Ar— 1, we consider the BAMP polynomials Pq'^ with L = m (mod N). Set 



b{t) = (Eilo^ *')"" in (4), then Fij{s) = P^'{s) with L = Nl + m and j = aL,Q. Hence (12) follows 



from (5), so wc obtain 

Suppose that Ak{s), 1 < k< N, are polynomials such that the relation (5) holds for Pq 



for all L and Q. By (11) and the definition of the matrix R in (2), one concludes that R^^^ + 



N\k 



0, hence Ak{s)'s are defined by the relation (3). As the matrix R 



N 



responds to the multiplication of {^j^)'^ on C[t]{= T,fjo't^C[s]), and (- 
J2qZi Fq^\s), one has the matrix-form expression of R^: 



-1\N 



t-i 



cor- 



R 



N 



( p(^) 
M) 

p(iV) 



,(iV) 



sP{ 



(N) 



AN) 



sP 



(TV) 
2 

(TV) 



AN) 



d{N) 
N-1 



sP 



(N) \ 



N-1 



piN) 



whose trace gives Ai{s) = —NPq^\s). By detii = (1 
= (— l)^(s — 1)^(^~^). It is easy to see that Pq^\s) is a monic polynomial of degree {N — 1) with 
Pq^\o) = 1, and Pq^\s) for Q 7^ are polynomials of degree {N — 2). Therefore the {i,j)th. entry 
of R^ is a polynomial of degree A^ — 1 or A — 2, according to i < j or i > j respectively, which 
implies degAk{s) < k{N - 1). Setting s = in (3), one obtains {x - 1)^ = + X;f=i Ak{0)x^ 



\N-1 



one obtains Ai\[{s) = (— 1) det R 



TV 
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hence ^a;(0) = (-l)'^(^). Substituting s by s ^ in (12), then multiplying the factor s''i+Jv2,Q^ by 
(i) we obtain the following relation for L and Q' with L + Q + Q' = (mod N) : 

p(^^')(s) + . . . + s^(^-i)Afe(i)p(^(^-'=)^)(s) + . . . + s^(^-i)A^(^)p(^)(s) = . 

As the above relation of BAMP polynomials holds for all L and Q', by the characterization of 
Ak{s) in (a), one obtains s''^^"^^ Ak{-) = Ak{s) for 1 < k < N. Therefore Ak{s) is a reciprocal 
polynomial of degree k{N — 1) with the leading term equal to ^^(0) = (— 1)'^('^). As the roots 
of a reciprocal polynomial consists of pairs {si,s^^} for Sj 7^ ±1, together with certain possible 
multiple-roots at s = ±1, the conclusion for Ak{s) , A; > 2, in {in) follows from (7) and the reciprocal 
property of Ak{s). □ 

Remark Since the coefficients of Pq'\s) are all positive integers, its real roots must be negative. 

By the relation between P^\s) and the eigenvalues of the Hermitian operator (1) (which will be 

explained in the next section), all the roots of Pq'\s) must be real numbers through the Onsager's 
algebra representation theory. However, a mathematical argument for this real-root property of 



Pq"^ (s) purely from the polynomial-algebra viewpoint has not yet be found. 



□ 



By (6), we list the recurrence relation of BAMP polynomials for iV = 2, 3,4 for later use: 
N = 2, P(^+^) - 2is + 1)P(^-^^) + is- l)^Py) = 0, 

N = 3, pJ-^+^) - 3(s2 + 7s + l)Pf + 3(s - l)4p^^+3) -{s- \fP^^^ = 0, 



iV = 4, P^^^^^ - 4(s + l)(s2 + 30s + 1)P<^^+^^^ + 2(3s2 - 62s + 3)(s - ifP^Q^^^ ^"^^^ 
-4(s + l)(s - 1)«P(^+^) + (s - l)i^Pf ) = . 

For a general iV, there is a general expression of Ak{s) as follows. 

Proposition 1 Denote Xg = J2q=o 00'^^ for < £ < N — 1. Then the polynomials Ak{s) for 
1 < k < N in (12) are expressed by 

Ms) = {-!)' E ^^'^•••^4- 

^i«?2<...<4 

Proof. With R in (2) and D := dia[s, , s^], D~^RD is the N x N Toeplitz (cychc) matrix 

with s^izi.)/^ as the {i,j)th entry. Hence is the eigenvalue of ( s^izZ)/^ J 

0<i,j<N-l \ J i,j 

with the eigenvector {l,uj^ , . . . ,uj^^~^^^Y for < ^ < N — 1, and D~^R^D is conjugate to 
the diagonal matrix dia[Ao^, • • • , A^_^]. Since R^ and D'^R^D have the same characteristic 
polynomials, the expression of Ak{s) follows immediately. □ 

Remark. In the A^-expression, s^^^ is the variable t. Under the transformation t ^ cut, Xi is 
changed to Af_i with A_i := Ajv-i- Therefore the A^-expressions in the above proposition do give 
the s-polynomials, which are equal to Ak{s). □ 



4 Differential Equation of BAMP Polynomials 

In this section, we discuss the differential equations satisfied by BAMP polynomials Pq"^ ■ For each 
L, we denote by V^^^s) the A'"- vector with Pq^ as the jth component for < j := o^^q < N — 1, 
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i.e., V^^\s) is the A/"-vector, 



Then with the matrix R in (2), we have 



Af-l 



(14) 



(15) 



j=0 



By differentiating t-^ the first relation of (15), then multiplying it by J2jLo^ we obtain 

(l,t,---,t^-i)-iVs(s-l)p(-^)'(s) 
= (l,t, • • • ,t^-') • (i^(« - l)(Ef=o'ii^) -{s- l)dia[0, 1, • • • , (iV - l)])p(^)(s), 



(16) 



where we use the superscript prime ' to denote the differentiation with respect to the variable s for 
the rest of this section. By 

N-\ N-l 
j=0 j=l 



sL 



Jn-, 

the relation (16) is equivalent to the following differential equation of P^^\s), 

Ns{s - i):pw = h:pW 



, 1 < J < iv - 1 , 



(17) 



where B = L{N— l)s — LY,j^i Bj — (s — l)dia[0, 1, ■ ■ ■ , (N — l)], which has the following expression: 



B 



( 5q —Ls ■ ■ ■ —Ls \ 
-L Si '■■ : 



V 



-L 



•. -Ls 

-L Sn-1 J 



, Sj := L{N - l)s - j{s - 1) . 



By successive differentiations of ^ on (17), one can express A/"*s*(s — 1)^^^^^ in terms of P^^\ 
In particular, for k = 2,3, we have 

N^s'^is - ifV^L)" = (b"^ - N{2s - 1)B + Ns{s - l)B'^p(^^ 

Ar3s3(s _ i)3-p(L)/// ^ Ib^-3N{2s-1)B^ + 2N^{Ss^-Ss + 1)B (18) 

+Ns{s - 1){2B'B + BB') - 2N^{2s - l)s{s - 1)B'^P^^^ . 

Using the relations between higher-order derivatives of V^^\ the following statement is expected 

to be true. 

Conjecture 1. There exist diagonal matrices {0 < k < N — 1) with entries (depending on 
L) in C[s], such that 



(19) 



k=l 



By which Pq"^ satisfies a Nth order differential equation with the regular singular points at s = 0, 1, 

and the expression depends only on L + Q (mod N). □ 

We are going to demonstrate the solution of (19) for N = 2,3. For A'" = 2, by (17) and (18), 
one has 



2(s - l)sP(^)' = 

4(s - ifs^p(L)" 



p{L) 



Ls —Ls 
-L Ls-{s~ 1) 

s(s + 1)L2-2s2l' -2s2l2 + s(3s - 1)L \ p{L) 

-2sL2 + (5s - 3)L s(s + 1)L2-2s(2s-1)L + 3(s-1)2 J 



Then one obtains the solution of (19) for N = 2: 

4s(s - l)p(^)" + 2Dip(^)' + DoV^^^ = 0, 
with Di = -dia [/o, /i] , Dq = dia [go,gi] for 

/o = (2L - 3)s + 1, /i = (2L-5)s + 3; 5o = i^(i^-l), ffi = (L - 1)(L - 2) . 
The differential equation of Pq"^ is given by the hypergeometric differential equations, 
L..en : 4s{s - l)P<f - 2f,_^^ (.)pf > + g,_^^ (,s)pf ^ = 0, 

with the polynomial solutions: Pq'\s) ^ F{^^^^—^,^y^,Q + ^;s) for even L, and Pq'\s) 

P(l - i - I - Q; s) for odd L. 
For N = 3, (17) and (18) become 



3s(s - 1)P(^)' = P(^); 9s2(s - l)27?(i')" = (^-^ _ 3(2s - 1)B + 3s(s - 1)S'^ p(^); 



{b^ - 9(2s - 1)^2 + 18(3s2 - 3s + l)B + 3s(s - 1)(2B'B + HB') 
■18(s - l)s(2s - 1)^') P(^) , 



with 



/ 2Ls 



B = 



-Ls 



-Ls 

L 2Ls - (s - 1) -Ls 
L -L 2Ls - 2(s - 1) / 



( 2L -L -L \ 

2L-1 -L 
2L - 2 y 



The solution of (19) for AT = 3: 

27s2(s - 1)27?W" + 9s(s - \)D2V^^^" + SD^r^^y + Do:P^-^^ = , 
is given by D2 = -3dia [/o, /i, /2], = dia [50, 51, 52], Do = -{L- l)dia [/iq, /ii, with 

/o = (2L - 4)s + 2, £/o = 3L2s(4s - 1) - 3Ls(10s - 7) + 2(s - l)(10s - 1), 

/io = i^(M8s + l)-4(s-l)^, 
fl = (2L - 6)s + 4, gi = 3L'^s{As - 1) - 3Ls(18s - 15) + 2(s - l)(31s - 10), 

hi = {L- 2) (^L{8s + 1) - 12(s - 1)^ , 

/2 = (2L - 5)s + 3, 52 = 3L2s(4s - 1) - 3Ls(14s - 11) + 2(s - l)(19s - 4), 
h2 = L'^{8s + 1) - L(16s - 7) + 6(s - 1). 
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Then one obtains the third-order differential equation for Pq for L = (mod 3): 

27s\s - l)2p(^)'" - 27s{s - l)fQPi}^" + ^gqP^^^' - (L - l)hQP^Q^^ = , 
The above differential equation holds also for P^^ for a general L with L + Q' = Q (mod 3). 

5 Chebyshev-type Polynomials nQ^(c) 

For the eigenvalue problem of (1), one employs the theory of Onsager's algebra representation 
through the zeros Cj's of some polynomials in the variable c, which relates to the variable s in the 
previous two sections by the transformation, c = jz^, equivalently, s = The change of variables 
provides the one-to-one correspondence between ,s G (— oo,0) and c G (—1, !)• Accordingly, for the 
description of eigenvalues of (1), one may convert a s-polynomial F{s) to a c-polynomial n(c) via 
the relation 

n(c) := (c+l)deg^F(^). 
By which the c-polynomial corresponding to the BAMP polynomial Pq"^ (s) will be denoted by 

n(3^)(c) = (c+l)^-.'?p(^)(^). (20) 

Then the eigenvalue of (1) (of size L) for the lowest state eigenvector in the Q-scctor is given by 
E = 2Qk' + {NbL,Q -{N-l)L){l + k')-N Y^jt? T^Jl + 2k'cj + k'"^ , where cj arc zeros of n[^^ (c) 

[3]^. We shall call IIq'\c) the (generalized) Chebyshev-type polynomials for the reason being clear 
later on in the discussion of the case N = 2. 

The reciprocal relation of s-polynomials Pq'\s) is translated to the *-relation: IIq'\c) 

IIq'^*{c). Here p*{c) is the polynomial defined by p*{c) := (— 1)*^^^ ^p{—c) for p{c) G C[c]. Hence a 
*-symmetric polynomial p{c) ( i.e., p{c) = p*{c)) simply means an even or odd function according 
to the parity of the degree of p(c). Then the relation , (pq)* = p*q*, holds. However, the equality, 
{p + q)* = p* + q*, is valid only for p, q with the same degree. 

Lemma 1 nQ''*(c) is a polynomial of degree b^^Q with the leading coefficient equal to N^~^ . And 
nj^^* (c) = H^q) (c) if and only if L + Q + Q' = (mod N) . 

Proof. By the positivity of all coefficients of Pq'\s), nQ'-*(c) is a c-polynomial of degree the same 
as the s-polynomial Pq'\s), which is equal to (Hence all the roots of H^^^ are real, confined 

in the open interval (—1, !)■) The leading coefficient of Hq^^(c) is given by P^^\l) = N^~^. The 
*-symmetry relation of Hg'^s follows from Theorem 1 {i). □ 

As a corollary of Lemma 1, one obtains the following results: 

Lemma 2 (i) When N is odd, for each L there is exactly one ^-symmetric Hq'^ with Q = 1^]^ 
(mod N). 

(a) When N is even, Hg'^ = Hg'^* holds only when L is even, in which case there are exactly 
two * -symmetric Hq'^ with Q given by the relation Q + ^ = (mod ^)- 

^The "coaOj" in formula (31) of [3] is equal to "— Cj" here due to the convention, c = yrf , used in this aricle, also 
see formula (33) of [16]. 
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□ 

Denote the c-polynomial associated to Ak{s) in Theorem 1 (m) by 

a,(c) := (c + l)M^-i)^fc(^), 1 < A: < iV . (21) 

By bi^kN^Q = bi^Q + k{N — 1), the substitution s = in (12), together with the multiphcation 

of (c + l)''i+Jv2,Q^ gives rise to the recurrence relation of IIq'\c) for a fixed Q. Indeed by Theorem 
1 (a) and {in), we have the following result: 

Theorem 2 With the polynomials ak{c) in (21), we have the following recurrence relation o/Hg'^s 
for the same Q, 

ng^+^^c) + . . . + a,(c)ng^+(^-^)^)(c) + . . . + a^m^^\c) = . (22) 

The polynomials ak{c) are *- symmetric, andai{c) = —NIIq^\c), aAr(c) = {—l)^2^^^~^\ afc(ibl) = 
(-l)'=(±2)'=(^-i)(^). For k>2, the degree of ak{c) is less than k{N - 1) with k{N - 1) - dega^ 
being a positive even integer. 

□ 

For AT = 2, 3, 4, the recurrence relation (22) takes the following explicit forms: 

N = 2, n(f+^)-4<+^)+<)=o 

N = 3, U^^^'^ - 3(9c2 - mii^'^ + ^8U^i^'^ - 64nf ) = , 

N = A, ng^+^«) - 32c(8c2 - 7)nf - 128(140^ - 17)ng^+«) - 2048cnf ^ ^ 
+4096n = 

In above discussions, the recurrence formula (22) follows from that of BAMP polynomials. 
However, it is more convenient to derive (22) directly from relations between IIq'){c)'s as follows. 
Parallel to (14), we define the following AT- vector 11^^^ and N x N matrix 3fi: 

/c+1 c-1 ••• c-l\ 



n(^) = (n^(c),n^(c), . . . ,0^3^(0))* , ^ 



c+l c+1 ■•• : 
: ■•. ■-. c-1 

\ c + 1 ••■ c + 1 c + 1 y 



Denote 

C{L) ^ a('^+i)sRA(-^)-1(c + 1)-^ (24) 

where A^^) := dia [(c + 1)^^'^, (c + 1)^^--^, . . . , (c + ifL,-L-N+ii^_ Foi < i,j < N - 1, we 
have b^+i-L-i-i > bL-L-j-, and the equality holds only for certain (i,j)'s with i > j. Hence 
the entries of C^^^ are all c-polynomials. Furthermore by bL+N,Q = ^L,Q + — 1), one has 
A(L+7V) ^ (c+ 1)^^-1 A(^). This implies C^^^ = C(-^+^), equivalently, C^-^) = C(^). By (15), one has 
the following relation, 

n{L+i)^^(L)n{L) (=c(^)nW), (25) 

hence n(^+^) = C^^+^-i) • ■ ■ C^L+i)^(L)^iL) ^ (c + l)-A'A(^+^)K^A(^)-in(^). Then we obtain 

YliL+N) ^ ^ -L^-l^(L)g^iV^(L)-ln(i) . (26) 
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Note that by bNk,Q = {N - l)k for all Q, U^^+^) = {c+ ly^^^U^^^ for L = 1 (mod N). The 
characteristic polynomials of the linear map in the right hand side of (26) is given by det(a: — (c + 
l)~^3f?^) = + J2k=i (^kic)x^~^ , which provides the coefficients in the recurrence relation (22) 



of n^^^'s. 



We now consider the case N = 2. By Lemma 1, Ilg 

(2 

Q 

t 



(2k)* 



degn[,^*^^ = degnf + 1 = degn^^+^^ = k. We have 



^(2fc) -p.(2A;+l)* _ ^(2fc+l) 



and 



jj(2fe) 



n 



(2fc) 



c),n, 



{2k). 



and C(^)'s in (24) arc described by 



, n(2'=+i) = (nf+^)(c),nf+^)(c) 

' c+1 c - 1 



1 c- 1 
1 c+1 



I . Then one has 



1 1 

n(2fe+l) = C^O)jl('^k) ^ jj(2fe+2) ^ ^(l)jj(2fc+l) . jj(2fe+2) ^ ^(1)^(0) jj(2fc) ^ jj(2fe+3) ^ ^(0)^(1) jj(2fc+l) 

1 \ , „ / C C+1 



The matrices, C(i)C(°) =2 . ^ ^ 

\ - 1 C 



C 

c- 1 



and CWC(I) = 2 

equation, x"^ — 4cx + 4 = 0. Hence one has the following recurrence relation of ng''"s 



, both satisfy the 



n(^+^)-4cn(.^+^)+4ng^) = o. 



which is the same as in (23) for N = 2. By normalizing Hq''' by a 2-power factor, one obtains the 
Chebyshev-relation: 

Proposition 2 For e = 0, 1, denote by Fj^ := 2^'^nQ^^'^'' for k G Z>o. The polynomials Fk{c) 
satisfy the Chebyshev recurrence relations: F^+i — 2cFk + Fk-i = 0, with Fq,Fi given by 



e 


Q 


Fo 


Fi 








1 


c 





1 





1 


1 





1 


2c + 1 


1 


1 


1 


2c- 1 



As a consequence, = 2 ^IIq^'^'', 1/^ = 2 '^n^^'^^ are the Chebyshev polynomials of first and second 
kind respectively^, and 2-'=nJ^''+^^ = Uk+i + Uk, = Uk+i - Uk- 

□ 

It is known that Chebyshev polynomials form a system of orthogonal polynomials satisfying a 
second-order differential equation. Indeed by using the relations between n^*^"'"^ , n^*^"*"^ and 11^*^, 
one obtains n{,^'=+^^ = 2cnJ^''^ + 2(c2 - l)nf and = 211^^^^ + 2cnf Then by induction 

argument, one can show 



{2k) 



dc 



(i-c^ 



dc 



-mf)+cnf) 



'Here we use the standard conventions: Tk{x) = cos(fcarccosa:); Uk(x) = sin(fcarccosx)/Vl — x'^- 
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which are equivalent to V1^£u^q''^ = kVl^uf''\ Vl^£{Vl^uf''^) 



-kU, 



(2fe) 



Both IIq ' and Vl — c^H^ are solutions of the equation, (Vl — c^^)^/ = 'i^^/; then follow the 
differential equations of IIq^, 



dc 



, (1 



^2TT(2fe) 



ac 



{2k) 



By which, nQ^'^^(c) for /c G N form a system of orthogonal polynomials with weight (1 — c^) 2 ; and 
the same for nj^^'^''(c)'s. 

6 Differential Equation of Chebyshev-type Polynomials 

Along the path of discussions in the case N = 2, one tends to find the differential equations of 
Chebyshev-type polynomials IIq''' for iV > 3, especially for those *-symmetric ones described in 
Lemma 2. By differentiating (20), one obtains the relation: 



dU 



Q 



- 1)^^ = bL,Q{c - - (C + l)l+^-.Q.(. - 1) 



dP, 



(L) 



Q 



dc 



ds 



Using (17), one reaches the system of differential equations of Hq'^'s: 



dc 



(27) 



where V = (c/ij)o<i,j<iv-i is the N x N matrix with entries given by 



-(c- 



,1+bL, 



hi -L-i ) - I7 for i 



^(c-l)(c + l)''-'^ 



-L-J 



for i > j, 
for i < j. 



Note that one has 



0, 
1, 



i < 1 



L<j,i, or j,i<)_ 
L<j, 



-1, j < 1 - L < i 



which imply entries of V are all c-polynomials. In particular, (27) for A'^ = 3 becomes 



(c^ 



(c^ 



dc 



2 

(3fe+l) 



V no 



(3fe+l) 



(3fc+2) 



^(3fe+2) 
jj(3fe+2) 



/ 



c'-l 

1 l-3c 
^ 3fc 
1 C+ 1 



'-1 

- 1 
l-3c 



/ nf )(c) \ 



3fe 



n 



{3k) 
2 I 

r(3fe) 



V nr^(c) J 



3k + 1 



3k + 2 



/ -2(c-l) -, -, 

/ ' c — 1 c — 1 



3fe+l 
C + 1 

V c+1 

c+1 



3fc+2 
C+1 
1 



-2c 
3fc+l 

C+1 

C- 1 
-4c+2 
3fc+2 
1 



C- 1 

-2(c+l) 
3fc+l 

c2-l 

c2-l 
-c-3 
3fc+2 



/ nf +^)(c) \ 

nf+^)(c) 
V nf -^^^(c) J 

nf +2)(c) \ 



n 



(3fc+2) I 
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Due to the complicated expressions as indicated in above, it seems quite difHcult to obtain a general 
form for the differential equations of Hq'^ for an arbitrary N. However for L = (mod N), the 
relation (27) becomes 



n(^) 



n 



V 



L 

N 



I 



( 
1 



V 1 



-2-Ar(c-l) 
L 

C+1 



C+1 



c- 1 

-4-Af(c-l) 
L 



C+1 



C^-l 
C- 1 



c- 1 

-2{N-\)-N(c-\) 
L 



\ 



( n(^) \ 



n 



i.L) 
N-l 



The first component in the above system gives the following simple relation: 



dc 



J N-l 

-^ng^) for L^O (modiV) 

Q=l 



(28) 



By using (27), one expects the following statement to be true for all N as in the case N = 2: 

Conjecture 2. The polynomial ^q{c) satisfies a ATth-order differential equation with regular 
singular points at c = ±1. □ 

We are going to justify the above conjecture and derive the differential equations of *-symmetric 
polynomials IIq^^ (i.e., IIq'^''''^^) for N = 3. For the rest of this section, we consider only the case 

N = 3 where the degree bL,Q of IIq'\c) are given by 

2fc = 63^,0 = ^3fe,l + l = ^3fe,2 + l = hk+lfi = ^3fe+l,l = &3fe+l,2 = &3fe+2,0 — 1 = &3fe+2,l — 1 = ^3fe+2,2- 

By (23), one has the four-term recurrence relation of Hq^^s : 



64n, 



(3fc+e) 



, 



(29) 



with IIq^^'^'' for A; = 0, 1, 2 given by 



e 


Q 






Q 








1 


9c^ - 5 


243c^ - 270c^ + 43 





1 

2 





9c ±3 


243c3 ± 81c2 - 135c + 21 


1 



2 


1 


27c2 ± 18c - 5 


729c^ ± 486c3 - 540c2 + 270c + 43 


1 


1 


1 


27c2 - 11 


729c^ - 702c2 + 85 


2 




1 


3c ± 1 


81c3±27c2-57c=Fll 


2187c5 ± 729c^ - 27540^ + 702c2 + 711c ± 85 


2 


2 


3 


81c2 - 21 


2187c^ - 1782c2 + 171 



By Lemma 2 (i), 

„{3k+Q)* _ Tr(3fc+(5) T-r(3fc)* _ ^{3k) „(3fc+l)* _ Tr(3fc+1) „(3fc+2)* _ Tr(3fc+2) 
~ ' ^^2 — ^^1 5 ^^2 — ^^0 ' ^^1 — ^^0 

Hence one can express the relation (25) in terms of Hq*^"'"'^^ and H^^'^^ , Hq^*^"*"^-* , Hq^'^"'"^^ : 



Hi 



(3fe+l) 




= Hf) + (c + i)(nl' 



(3fc) 



H 



(3fc)*N 



jj(3fe+i) ^ nf ) + (c - i)Hf ) + (c + i)h(=^^> 



jj(3fc+2) ^ _ ^^j^(3.+l) ^ ^ l)(nf +1)* + Hf +^)), 



(3fc+2) 
2 

(3fc+3) 




n, 



(3A:+1) 







(3A;+1) 



+ n, 



(3A:+1)* 




(30) 



i)nf+2) + (c + i)nf+2) 



jj(3fc+3) ^ +2) ^ +2)* + + l)nf . 



+ (c2 - i)n 



(3fc+2) 
2 ' 
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By (26) with (iV,L) = (3,3A;), one has 



(3A:+3) 



n: 



(3fe+3) 



9c^-5 3(c2 - l)(3c- 1) 3(c2 - l)(3c + 1) 
3(3c + l) 9(P-5 3(c+l)(3c-l) 

3(3c - 1) 3(c - 1) (3c - 1) 9c2 - 5 



(.3fc) 



Tim 



For convenience of notations, we denote u'^''\Hk := U.^i''\ One has 

Pfe+i = (9c2 - 5)Pjk + 3(c2 - 1) (^(3c - l)Hk + (3c + 1)H*^ , 

Hk+i = 3(3c + l)Pk + (9c2 - 5)Hk + 3(c + l)(3c - l)Hl 
iJ^+i = 3(3c - l)Pk + 3(c - l)(3c + l)Hk + (9c2 - 5)H^ . 

Lemma 3 Denote by = — ■ Then the following relations hold, 



(31) 



k{Hk + HI) = ^, M3|±l)i^, = _(c2 _ 1);^ + (fc _ i)c^ + 2k^Pk, 

2fc^(fc+l)c p 



1 



— 3 dT - i'^ + ^ic-^ - Ufc - ^; - gpniyj-ar ^-i 



(32) 



Proof. The first relation in (32) follows from (28) for A'' = 3. Wc now show the other two relations 
by induction on k. The relations are easily verified for = 1. Assume (32) holds for k. By (31) 
(32), we have 

Pfe+i = (9c2 - 5)Pfc + f c(c2 - 1)^ - 3(c2 - l)Ek; 

%i ={k + l){18cPk + l{l8c^-8)^-6cEk); 

Ek+i =6Pk + 6c{Hk + H^)-2Ek = 6Pk + - 2Ek. 

By differentiating the 2nd and 3rd relations in above, one obtains: 

_c!^<i!^ = (9(2A: + l)c2 -8k- 9)Pk + ic(9(2/c + l)c2 - 16A; - 13)^ 

-i(9(2A; + l)c2 - 12A; - 13)^;^; 
(c^ - l)%i = -UkcPk + |(-9c2 + 6^)^^ + AkcEk . 

By which follow the second and third identities of (32) for A; + 1. □ 



Theorem 3 For L = Q (mod 3), the polynomial Hq'^ satisfies the following 3rd order differential 
equation: 



27(c2 - 1) (^(c^ - 1)^^ + 2(Q + 2)c 



dU 



dc'^ 



9(L -2Q- 2)(L + 2Q + Z){c^ - 1) 



-12(3g2 + 3Q + 2)j - 2{L - Q){L + 2Q)(L + 2Q + S)cU'^q^ = . (33) 

Proof. In this proof, Pk,Hk,Ek denote the c-polynomials in Lemma 3, and the c-derivatives of a 
function /(c) will be denoted by /',/",••• etc. Among the relations in (32), the differentiation 

of the second one is the same as the left hand side of the third one, by which one obtains the 
constraint of P^, i.e. IIq^'^^ which is the differential equation (33) for Q = 0. 
By (30) and (32), we have 



(3fc+l) 



(3A;+1)/ 



Pk + tPL - Ek 



Pfc + c{Hk + H^) — Ek — ^ k -r k 
2kcPk + {2{c'-l)-i^)P;^ + ^Ek, 



(34) 
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which impUes 

2cnf ^+1) + 3(c2 - l)nf *^+^)' = 2{3k + l){cPk + ^^Pi) , 

m - 'fEk = -{c' - l)nf + 2k{cP, + ^^Pl) (35) 

- (c"-l) rr(3fc+lV I 2kc yfi^k+l) 
~ 3fc+l ^^1 3fc+l^^l 

Differentiating the 2nd relation in (34), then using the 2nd one in (32), we have 



(c2 - l)np'=+')" + 2cnf + (^P^ - f Efe)' = 2kPk + 2{k + 2)cP^ + 2(c2 - 1)P^' 

= 2A;(2A; + l)Pfe + 2(2A; + l)cP^ - 2^(3^+1) ^^ . 



By the last equality of (35), we obtain 

2(21: I IIP.. I ^(^^ + ^)" P^ 2(3fc + l) 3(c^-l) (3fe+i). 8c (3.+1)/ 2 (3,+i) 
2(2fc + l)Pfe+ ^ P, ^ ^,____n, +SkTl^' +3fcTT"^ 

Then from the first relation of (34), one arrives the following relation: 

2(3fc + 2) 2(3fc + 2)c , _ 3(c2-l) ^(3fc+i)// 8c ^(3fc+i)/ _ 2(9P + 6fc - 2) ^(3^+1) 
3 3fc 3fc + l ^ +3fc + l 1 3(3)t + l) 1 

By which and the first relation in (35), one solves Pfe,Pfe in terms of nf +\ ilf^^^' and Uf^^^": 

-2{3k + l)i3k + 2)Pk = 9(c2 - l)2nf '^+^)" + 9(2 - A:)c(c2 - 1)0^+'^' 

+2(9A;2 + Qk-2- 9kik + l)c2)nf , 
2(3k+i)(3k+2) p, ^ 3^^^2 _ i)nf + 3(2 -k)cP + {3k + 2)nf - 6k{k + l)cnf . 

By comparing the derivative of the first equality with the second relation in above, one arrives the 
differential equation (33) for Q = 1, i.e. 11^*^"^^^ 
By (30) and (34), we have 



(36) 



(c2 - l)nf = 6kcPk + (6c2 - M)P^ - 2kcEk , 
which implies 

^^=-2(3^ + 2) 4(3^ + 2)^"^ ^^"2 • ^^^^ 

By differentiating (37), and using the second relation of (32), we obtain the following identity from 
expressions of (c^ — 1)P^': 

2PP, + ik- i)cP^ - = -^(c^ - D^nf + f^c(c^ - i)nf +^)' 

+ 2(3fc+2) ~ ■'■^^2 

Then by (37) and the first relation of (36), we obtain 

12(A;+l)(3A:+2)Pfc = 9(c2-l)2n(^+^)"+36c(c2-l)n^^+^)'+2f (3fc+4)(6A:+l)-9A;(2A;+3)c2')n(^+^) . 



By comparing (37) with the derivative of the above relation, one obtains the differential equation 
(33) for g = 2. □ 
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The *-symmetric Hq'^"'"'^^ can be regarded as N=3 version of the Chebyshev polynomials ap- 
peared in the case N=2. The rest can be expressed in terms of 11^'^''"*^^ 's. Indeed, by Lemma 
2 [i), IIq'^ for L ^ Q (mod 3) are determined by the following *-symmetric or antisymmetric 

polynomials: nf ^ - nf )*,nf + nf +^)*,nf - nf +^)*, and nf ^ + nf )*,nf - 

jj(3fc+i)* jj(3fc+2) ^ nj3fc+2)* Lemma 3, we have the relations, 



jj(3fc) ]-[(3fc)* 
„(3fe+l) y^(3fc+l)* 
l-r(3fe+2) y,(3fc+2)* 

IIq -t- ilg 



(3A:+2) 



(3k) _ jj(3fc)* 



n(3fc+i) 



4dn^ 

k dc 



{3fc) 



n 



jj{3k+l) (3fc+l)* 

T-r(3fc+2) ^(3fc+2)* 
^^0 ~ ^^0 



2^(3^-1). 

2rfn^. 

k dc ' 

2n(3'=+i) . 



By which, one can express the functions in terms of n^'^'*"'^^'s. 

Remark. As n[J''+'^^ is an even c-polynomial of degree 2k, hence depending only on c^, one can 

{3k+Q)^^-^ 



express IIq''^'^\c) as a polynomial of the variable a 



TTQ\a) for some 



cr-polynomial -Kq^ of degree k, with all its roots in the open interval (0, 1). By (29), tTq^ satisfy 



the 4-terms recurrence relation: ttq"''^'* — 3(4 — 9a")7rQ + 487r, 
A; = 0, 1, 2 given by 



(fc+i) 
Q 



647r, 



^■■^ = 0, with vrg'^ for 



Q 






_(2) 

Q 





1 


-9cr + 4 


243cj^ - 216cr + 16 


1 


1 


-27a + 16 


7290-2 - 756cr + 112 


2 


3 


-81C7 + 60 


21870-2 - 2592cr + 576 



Using the relations, £ = -2c A, |, = 4(1 - a)^ - 2^ and ^ = -4c (2(1 - a)^ - 3^) , 

(k) 

can convert the equation (33) into the following differential equation of ttq : 

(j2 {*) J (fc) 

4a2(l - a)^ - 2a((l + 2Q)a - 2Q)^ + ((3^^ + (2Q + l)k - + 3Q - 2)a 

: 4(3Q^-6Q+2) )^W^ = -(it - 2 + Q)(fc - 1 + Q)(A; + l)7r(J) . 



one 



-1 



□ 



Acknowledgments 

This paper is mainly based on the joint work with G.von Gehlen [14], except that a systematic 
account on mathematical structures of BAMP and Chcbyshcv-type polynomials for an arbitrary N 
is presented, and certain parts contain variations from the theme in [14] by adding more rigorous 
mathematical arguments here. This is an occasion to thank G.von Gehlen for the stimulating 
collaboration and correspondences. This work has been supported in part by NSC 93-21 15-M-OOl- 
013, Taiwan. 



References 



[1] G. Albertini, B. M. McCoy, and J. H. H. Perk, Eigenvalue spectrum of the superintegrable 
chiral Potts model, in Adv. Stud. Pure Math., 19, Kinokuniya Academic (1989) 1-55. 



17 



[2] R. J. Baxter, Exactly solved models in statistical mechanics, Academic Press, London, 1982. 

[3] R. J. Baxter, The superintegrable chiral Potts model, Phys. Letts. A 133 (1988) 185-189. 

[4] R. J. Baxter, Chiral Potts model with skewed boundary conditions, J. Stat. Phys. 73 (1993) 
461-495. 

[5] R. J. Baxter, Interfacial tension of the chiral Potts model, J. Phys. A: Math. Gen. 27 (1994) 
1837-1849. 

[6] T. S. Chihara, An Introduction to Orthogonal Polynomials, Gordon and Breach, Science Pub- 
lishers, 1978. 

[7] S. Dasmahapatra, R. Kedem and B. M. McCoy, Physics beyond quasi-particles: Spectrum 
and completeness of the 3 state superintegrable chiral Potts model, Nucl. Phys. B 396 (1993) 
506-540, hep-th/9204003. 

[8] E. Date and S. S. Roan, The structure of quotients of the Onsager algebra by closed ideals, J. 
Phys. A: Math. Gen. 33 (2000) 3275-3296, math.QA/9911018. 

[9] B. Davies, Onsager's algebra and superintegrability, J. Phys. A: Math. Gen. 23 (1990) 2245- 
2261. 

[10] B. Davies, Onsager's algebra and the Dolan-Grady condition in the non-self case, J. Math. 
Phys. 32 (1991) 2945-2950. 

[11] L. Dolan and M. Grady, Conserved charges from self-duality, Phys. Rev. D 25 (1982) 1587- 
1604. 

[12] G. von Gehlen, Onsager's algebra and partially orthogonal polynomials, Intern. J. Mod. Phys. 
B, 16 (2002) 2129-2136, hep-th/0201195. 

[13] G. von Gehlen and R. Rittenberg, Z„-symmetric quantum chains with infinite set of conserved 
charges and zero modes, Nucl. Phys. B 257 (1985) 351-370. 

[14] G. von Gehlen and S. S. Roan, The superintegrable chiral Potts quantum chain and generalized 
Chebyshev polynomials, in Integrable Structure of Exactly Solvable Two- Dimensional Models 
of Quantum Field Theory, NATO Science Series II 35, Kluwer Academic Publ. 2001, 155-172, 
hep-th/0104144. 

[15] L. Onsager, Crystal statistics. I. A two-dimensional model with an order-disorder transition, 
Phys. Rev. 65 (1944) 117-149. 

[16] S. S. Roan, Onsager's algebra, loop algebra and chiral Potts model. Preprint Max-Planck-Inst. 
fiir Math., Bonn, MPI 91-70, 1991. 

[17] G. Szego, Orthogonal Polynomials, Amer. Math. Soc. Colloq. Publ. 23, 3rd.ed,, New York, 
1967. 



18 



